Entanglement of two qubits mediated by one-dimensional plasmonic waveguides 
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We investigate qubit-qubit entanglement mediated by plasmons supported by one-dimensional 
waveguides. We explore both the situation of spontaneous formation of entanglement from an un- 
ent angled state and the emergence of driven steady- state entanglement under continuous pumping. 
In both cases, we show that large values for the concurrence are attainable for qubit-qubit distances 
larger than the operating wavelength by using plasmonic waveguides that are currently available. 
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As a direct consequence of the quantum superposition 
principle, a system composed of subsystems has states 
that cannot be factorized in products of states of its 
components. This non-separability, labeled as entangle- 
ment, is at the heart of quantum cryptography, quan- 
tum teleportation, or other two-qubit quantum opera- 
tions [TJ [2]. Exploited at first in systems like optics, 
atoms, or ions, entanglement is becoming more and more 
attainable in condensed matter physics. In particular, 
short distance entanglement is now available for spin or 
charge degrees of freedom in quantum dots (QDs), nan- 
otubes, or molecules [3H6]. However, for transmission 
of information at long distances, large separations be- 
tween the components are needed. For this purpose, cor- 
relation among the two qubits must be mediated by vir- 
tual bosons. Photons, either in the range of microwaves 
for coupling superconducting qubits [7] or in the visible 
range for QDs [8HTU], molecules, or NV centers in dia- 
mond [11 , are the usual candidates to play this role. 

Here we investigate a feasible proposal for long- 
distance entanglement of two qubits by using plasmons 
instead of photons. We consider the plasmon-polariton 
modes supported by one-dimensional (ID) plasmonic 
waveguides (PWs), see top panel of Fig.l. PWs have 
been studied during the last years as promising candi- 
dates to build up a new kind of photonic circuitry [T2] . 
The propagating plasmons associated with these struc- 
tures are characterized by both a subwavelength light 
confinement and long enough propagation lengths [13]. 
Coupling between quantum emitters and PWs has been 
also addressed p~4j [15]. These works show that the j3- 
factor, which measures the fraction of the emitted radi- 
ation that is captured by the propagating mode, can be 
close to 1 in realistic PWs. This is due to the subwave- 
length nature of the plasmon field in a 1D-PW. Very re- 
cently, these large /3-factors have been exploited to mod- 
ulate the energy transfer and superradiance phenomena 
appearing when two quantum emitters are placed at 1D- 
P Ws like channel or wedge structures [16] . In this Letter 
we show that PWs can also be used to obtain a large de- 



gree of entanglement in two qubits separated by distances 
larger than the operating wavelength. 
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FIG. 1. (Color online) (a) Two qubits interacting with a 
plasmonic waveguide, in this case a channel waveguide, (b) 
Scheme of levels, couplings, and decays in the particular case 
where uj\ — UJ2 — cjo and 711 — 722 = 7. 

The dynamics of the density matrix p for two qubits is 
described, after tracing out over the degrees of freedom 
of the plasmons, by a master equation [T71 [18] 

d tP =^H}+ ^faipai-vfap-pafa), (1) 

where cr},^ are the raising and lowering operators for 
each qubit. The ingredients of Eq. (1) are determined 
by the classical Green's function describing the electro- 
magnetic interaction between two dipole moments, \±\ 
and /i2, placed at locations 1*1 and r2, G(o;, ri, V2). For 
two qubits with the same characteristic frequency c^o, the 
hamiltonian can be written as 

H = hhJQ ^ a i a i + #12(^2 + 02<7i). (2) 
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The coherent part of Eq.(l) reflects the effective in- 
teraction between the qubits that is provided by the ex- 
change of virtual bosons [18]: 



912 



i r 



2 /m[^G(o; 7 ri 7 r2)/i2] 

C 2 (U - U ) 



duo, 



(3) 



whereas the rates of the non-coherent terms are given by 
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with i, j = 1,2 and 712 = 721- Equations (3) and (4) in- 
volve a point-dipole emitter approach, which is accurate 
enough for qubits such as atoms, small molecules or NV 
centers in diamond. For big molecules or QDs with sizes 
of a few tens of nanometers, a more realistic description 
of the quantum emitter is usually required [19] . 

When the plasmon supported by the PW is the dom- 
inant decay channel (i.e., large /?- factor), a very good 
approximation for the total Green's function can be 
obtained by only considering its plasmon contribution, 
G(cj,ri,r2) « G p i((jj,ri,r2) US]. In this way, analytical 
expressions for both gi2 and 712 can be easily derived: 



912 = \$e- d l^ sm(k pl d) 
712 = 7/^- d/(2L) cos(fc p id), 
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where k p i and L are the wavenumber and propagation 
length of the plasmon, respectively. These two magni- 
tudes, k p i and L, depend on the operating frequency,^- 
In deriving Eq.(5) we have assumed that the two qubits 
are equal and are placed at two equivalent positions along 
the waveguide, such that 711 = 722 = 7, and separated 
by a distance d. We define the modal wavelength of the 
plasmon, A p /, as X p i = 2ir/k p i. The crucial point of Eq.([5| 
is the 7r/2 phase shift between the coherent and incoher- 
ent parts of the coupling, which allows switching off one 
of the two contributions while maximizing the other by 
just choosing the inter-qubit distance. This opens the 
possibility of modulating the degree of entanglement. 

To test the feasibility of our proposal, we have carried 
out extensive numerical calculations on a particular PW, 
a V-groove milled on a silver film. In Fig. 2(a) we render 
the dispersion relation (energy versus wavenumber) of 
the propagating plasmon supported by a V-groove. This 
type of plasmon are usually called channel plasmon po- 
laritons (CPPs). The geometrical parameters are taken 
from realistic structures: the angle of the groove is 20° 
and its height is 140 nm, but similar results would be 
obtained for other sets of parameters. The evolution of 
the propagation length, L, of the CPP with the operating 
wavelength, A, is shown in the inset of Fig. 2(a). As ex- 
pected, L increases as A is enlarged. The dependence of 
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FIG. 2. (Color online) (a) Dispersion relation (red curve) 
of a CPP mode supported by a V-groove of angle 20° and 
height 140 nm. For comparison we show the dispersion of a 
surface plasmon mode of an infinite 2D silver surface (black 
dotted line). The inset displays the wavelength dependence 
of the propagation length, L. (b) /3-factor versus wavelength, 
A, and vertical distance, h, associated with the CPP mode. 



the /3-factor with both A and the vertical distance of the 
qubit(s) to the apex of the V-groove (see Fig. 1(a)) is dis- 
played in Fig. 2(b). Importantly, j3- factors larger than 
0.9 are attainable for a broad range of A's and within 
a large spatial region. Let us remark that the best f3- 
factors attained with dielectric waveguides are typically 
much lower than those obtained with PWs (for instance, 
/3-values not higher than 0.5 are reported in [16 for a 
GaAs fibre of 50 nm radius). 

Once we have introduced the ingredients of Eq. (1), 
now we discuss how this equation is solved. The most 
adequate basis to represent the dynamics of Eq.Q is 
the one depicted in Fig. 1(b): {|0) = |#i, #2) , |±) = 
^(|ei,# 2 ) ± l#i> e 2)),|3) = |ei,e 2 )}, where gi/a labels 
the ground/excited state of the z-qubit. Depending on 
both the sign and absolute value of 712, one of the states 
|±) can be practically decoupled from the dynamics of 
the rest of states. Once the density matrix p(t) is ob- 
tained by numerically solving Eq. ([T]), the entanglement 
of the two qubits is quantified by means of the concur- 
rence, C, defined as proposed by Wootters [20]. The 
two main ingredients controlling the dynamics of the two 
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FIG. 3. (Color online) Concurrence as a function of time when 
just one of the two qubits is initially excited. Continuous 
black line corresponds to two qubits entangled by means of a 
channel PW with /3 = 0.9, k p id = 2ty. Dashed red line is for 
two qubits in a cavity with a detuning such that it produces 
#12/7 = 5 and 712 = 0. Dotted blue line is for the ideal case 
of both PW and CQED (see text). The inset shows the time 
evolution of C and the populations of states |+) and |— ) for 
the non-ideal PW case. 

qubits (g\2 and 712) affect C in very different ways. The 
coherent coupling gi2 produces oscillations, whereas the 
cross-decay term, 712, produces a non-oscillatory contri- 
bution to C. These two effects are discussed below in two 
different situations. First, we analyze the case in which 
the system is initially prepared in a given unentangled 
state from which it decays spontaneously. In the second 
situation, the two qubits are continuously pumped by an 
external laser to reach a stationary state. 

In order to analyze the spontaneous formation of en- 
tanglement, one can initially prepare the system in the 
|ei, #2) state. In this case, the concurrence takes the form 

C{t) = y/[p ++ {t) - p__(t)]2 + 4/m[p + _(i)P, (6) 

where p±± are density matrix elements in the basis |±). 
The dynamics of C is shown in Fig. [3] for three different 
situations. With a solid black curve we render C(t) for 
the case of a channel PW where 712 and g±2 are given 
by Eq. (H), with fi = pe~ d ^ 2L ^ = 0.9 and k pt d = 2tt. 
The concurrence is characterized by a fast initial increase 
followed by a very slow decay. For this case, the coher- 
ent oscillations produced by g±2 are completely quenched 
(#12 = 0) and the cross-decay term (712) dominates. This 
dynamics can be easily understood from the time evolu- 
tion of the populations of the two entangled states |±) 
(see inset of Fig. 3). These states are equally popu- 
lated initially but the decay of the state |+) is very fast 
(7 + 712) while the decay of |— ) is very slow (7 — 712). 
The asymmetry between the two cascades is responsible 
for the long lifetime of C while the imbalance among the 



populations of states |±) determines C due to the first 
term in Eq.([6|. With a dotted blue line we represent 
the ideal case that corresponds to /3 = 1 (/? = 1 and 
L = 00) and h v \d — 2ir. In this case, the concurrence 
tends asymptotically to a steady state value of 0.5. 

It is worth comparing our PW-based entanglement 
with other schemes for achieving large entanglement that 
have been proposed before. In particular, embedding two 
qubits in a photonic cavity (CQED) [8 offers many possi- 
bilities for controlling the photon emission. As entangle- 
ment in CQED only relies in the coherent term, ujq must 
be tuned to the frequency of the cavity mode. In the 
case of perfect tuning [17 , the evolution of C with time 
is determined by 712 = # 2 /7, g being the qubit-cavity 
coupling. The time-dependent C becomes equal to that 
of an ideal PW (dotted blue curve in Fig. 3). However, 
in realistic implementations of CQED, one must work 
with a detuning A that is comparable to g. By using a 
Schrieffer-Wolf transformation one can obtain a master 
equation like Eq.Q with g±2 — g 2 j A and 712 ~ 0. In 
Fig. 3 we render a simulation to illustrate CQED-based 
entanglement (see red dashed curve) in which we have 
used #12/7 = 5 and 712 = as realistic values taken from 
experiments [9l[T0]. Although in this case C can be signif- 
icantly larger than the one obtained with realistic PWs, 
this only occurs in a very short time scale and presents 
very fast oscillations. Another advantage of the use of 
PWs is that they enable the emergence of long-distance 
qubit-qubit correlations, as we show below. 

A continuous pumping is required in order to have 
a stationary state with a high degree of entangle- 
ment. For sufficiently separated qubits, the stationary 
state can be modulated by acting independently and 
resonantly (cJi aS er = ojo on each qubit with a laser beam 
of Rabi frequency IV A new term, ^ Wti{a\ + <r^) must 
be included in the hamiltonian given by Eq.Q. When 
the system is initially prepared at state |0), one can ap- 
ply the continuous laser only on qubit 1. Inset (b) of 
Fig. [4] shows, for an almost ideal PW and very short dis- 
tances, the transient dynamics of the concurrence when 
d = Xpi/2 so that 712 is present while gi2 is quenched. 
At the beginning, clear oscillations are observed with the 
concurrence becoming larger than 0.5 for some time in- 
tervals, At ~ IO/7. 

More interesting is the discussion of the stationary con- 
currence, Coo, for different laser intensities and separa- 
tions between qubits. Here we consider realistic values for 
P and L taken from Fig. 2(b): f3 = 0.94 and L = 2 /mi, 
which correspond to a wavelength of around 640 nm and 
vertical distance h = 180 nm. Notice that for this h, 
the emitter-metal distance is about 47 nm, which is com- 
parable to the optimum distances found in other PW 
geometries like metallic wires or wedges [T6] , As shown 
in Fig. [4] (dotted blue line), when the system is prepared 
in the symmetric state |+) by a large spot laser which 
excites equally the two qubits (Qi = f^), Coo gets its 
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FIG. 4. (Color online) Panel (a) shows the steady state con- 
currence as a function of the separation between two equal 
qubits for f3 = 0.94, L = 2 /xm and three different laser con- 
figurations: Qi = 0.157 and Q2 = (solid black curve), 
^1 = ^2 = O.I7 (dotted blue line) and Qi — — Q2 — O.I7 
(dashed red line). Insets (b) and (c) correspond to the first 
type of pumping Qi = 0.157 and ^2 = 0. Panel (b) shows the 
concurrence time-evolution for a system with /3 = 0.99, while 
the /3-dependence of the steady-state concurrence evaluated 
at d = X p i is rendered in panel (c). 



maximum value for d close to an odd multiple of \ p i/2, a 
situation in which 712 is closer to —7 and the singlet state 
|+) is almost decoupled from the cascade decay involving 
the other three states (see Fig. [TJb)) . One can also pre- 
pare the state |— ) by pumping the two qubits with the 
same frequency and same intensity but phase-shifted by 
7r, which is equivalent to using ^2 = — ^1. As expected, 
the result is just shifted with respect to the previous case: 
Coo gets its maximum value for d being an even multiple 
of Xpi/2 (dashed red curve in Fig. 4). The X p i periodicity 
of these two previous cases changes to X p i / 2 when one of 
the two qubits is privileged with respect to the other by 
focusing the laser beam only on it (Qi ^> ^2)- Here we 
present the case ^1 = 0.157 while ^2 = (see solid black 
curve in Fig. 4). Notice that, in the three cases displayed 
in Fig. 4(a), large values for the concurrence at distances 
larger than the operating wavelength are attainable when 
using realistic values for both j3 and L. 

The robustness of our proposal with respect to changes 
in the /?- factor is analyzed in Fig. 4(c). As expected, the 
steady-state concurrence obtained at d = X p i is reduced 
when j3 is smaller than 1. However, this change is not 
very abrupt, as the concurrence is decreased by a factor 
of 2 when (3 is reduced from 1 to 0.8. Note that the 
scheme presented here for PWs could be also operative 
in other types of waveguides that display large /3-factors 
as, for example, photonic crystal waveguides [2Tj [22]. 

In conclusion, plasmon polaritons in realistic one- 
dimensional waveguides are excellent candidates to act as 
mediators for achieving large values of entanglement be- 
tween two distant qubits. We have shown that the large 
/3-factors associated with these waveguides and the fact 



that the coherent and incoherent parts of the two-qubit 
coupling driven by plasmons are phase-shifted allow to 
modulate populations and correlations between quantum 
states. We believe that our findings could also be useful 
for implementing the concept of dissipative engineering 
of states [2311251 - 

A.G.-T. and D.M.-C. contributed equally to this work. 
Work supported by the Spanish MICINN (MAT2008- 
01555, MAT2009-06609-C02, CSD2006-00019-QOIT and 
CSD2007-046- NanoLight.es) and CAM (S-2009/ESP- 
1503). A.G.-T. and D.M.-C acknowledge FPU grants 
(AP2008-00101 and AP2007-00891, respectively) from 
the Spanish Ministry of Education. 



[4] 
[5] 

[6] 
[7] 

[8 
[9 

[10; 

[11 

[12 
[13 

[14 
[15 

[16 



Corresponding author: carlos.teiedor@uam.es 
* Corresponding author: fj.garcia@uam.es 
M. A. Nielsen and I. A. Chuang, Quantum computing and 
quantum information (Cambridge, 2000). 
S. Haroche and J. M. Raimond, Exploring the quantum 
(Oxford, 2006). 

Y. Makhlin, G. Schon, and A. Shnirman, |Rev. MocT| 



Phys. 73, 357 (May 2001). 



R. Hanson, L. P. Kouwenhoven, J. R. Petta, S. Tarucha, 

and L. M. K. Vandersypen, Rev. Mod. Phys. 79, 4 (2007). 

J. R. Weber, W. F. Koehl, J. B. Varley, A. Janotti, B. B. 

Buckl ey, C. G. V. de W alle, and D. D. Awschalom, Proc. 

Natl Acad. Sci. USA| 107, 8513 (2010). 

S. D. Franceschi, L. Kouwenhoven, C. Schi'^nenberger, 

and W. Wernsdorfer, Nature nano. 5, 703 (2010). 

J. Majer, J. M. Chow, J. Gambetta, J. Koch, B. R. Joh- 

son, J. A. Schreier, L. Frunzio, D. I. Schuster, A. A. 

Houck, A. Wallraff, A. Blais, M. Devoret, S. M. Girvin, 

and R. J. Schoelkopf, [Nature 449, 443 (2007). 

A. Imamoglu, D. D. Awschalom, G. Burkard, D. P. Di- 

Vincenzo, D. Loss, M. Sherwin, and A. Small, Phys. Rev. 

Lett? 83, 4204 (Nov 1999). 

E. Gallardo, L. J. Martmez, A. K. Nowak, D. Sarkar, 
H. P. van der Meulen, J. M. Calleja, C. Tejedor, I. Prieto, 

D. Granados, A. G. Taboada, J. M. Garcia, and P. A. 
Postigo, [Phys. Rev. B| 81, 193301 (May 2010). 

A. Laucht, J. M. Villas-Boas, S. Stobbe, N. Hauke, 

F. Hofbauer, G. Bohm, P. Loda hl, M.-C . Amann, 
M. Kaniber, and J. J. Finley, [Phys. Rev. B| 82, 075305 
(Aug 2010). 

M. V. G. Dutt, L. Childress, L. Jiang, E. Togan, J. Maze, 
F. Jelezko, A. S. Zibrov, P. R. Hemmer, and M. D. Lukin, 
Science 316, 1312 (2007). 

T. W. Ebbesen, C. Genet, and S. I. Bozhevolnyi, Physics 
Today 5, 44 (2008). 

E. Moreno, S. Rodrigo, S. Bozhevolnyi, L. Martin- 
Moreno, and F. J. Garcia- Vidal, Phys. Rev. Lett. 100, 
023901 (2008). 

D. E. Chang, A. S. S0rensen, P. R. Hemmer, and M. D. 
Lukin, [Ph^s. Rev. Lett. 1 97, 053002 (Aug 2006). 
A. V. Akimov, A. Mukherjee, C. L. Yu, D. E. Chang, 
A. S. Zibrov, P. R. Hemmer, H. Park, and M. D. Lukin, 
Nature 450 (November 2007). 

D. Martin-Cano, L. Martin-Moreno, F. J. Garcia- Vidal, 
and E. Moreno, [Nano Letters|10, 3129 (2010). 



5 



[17] Z. Ficek and R. Tanas, [Physics Repor ts 372, 369 (2002), 
ISSN 0370-1573. 



[18] D. Dzsotjan, A. S. S0rensen, and M. Fleischhauer, Phys. 

Rev. B| 82, 075427 (Aug 2010). 
[19] M. Andersen, S. Stobbe, A. S. Sorensen, and P. Lodahl, 

Nature Physics (2010). 

[20] W. K. Wootters, [Phys. Rev. LetT] 8Q, 2245 ( Mar 1 998). 
[21] G. Le camp, P. Lalanne, and J. P. Hugonin, |Phys. Rev.| 
I Lett. 99. 023902 (Jul 2007). 

[22] T. Lund-Hansen, S. Stobbe, B. Julsgaard, 



H. Thyrrestrup, T. Siinner, M. Kamp, A. Forchel, 
and P. Lodahl, Phys. Rev. Lett. 101, 113903 (Sep 2008). 
[23] F. Verstraete, M. M. Wolf, and J. I. Cirac, |Nature Physics 
5, 633 (2009). 

[24] H. Krauter, C. A. Muschik, K. Jensen, W. Wasilewski, 
J. M. Petersen, J. I. Cirac, and E. S. Polzik(2010), |http: 

| //arxiv . org/abs/1006 . 4344| 

[25] A. F. Alharbi and Z. Ficek, |Phys. Rev.~A| 82, 054103 
(Nov 2010). 



